Critical Conductance of a Mesoscopic System: Interplay of the Spectral and 
Eigenfunction Correlations at the Metal-Insulator Transition 
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We study the system-size dependence of the averaged critical conductance g(L) at the Anderson 
transition. We have: (i) related the correction 5g(L) = g(oo) — g(L) oc L~ y to the spectral correla- 
tions; (ii) expressed Sg(L) in terms of the quantum return probability; (iii) argued that y = rj - the 
critical exponent of eigenfunction correlations. Experimental implications are discussed. 



PACS numbers: 73.40.Hm, 73.23.-b, 72.15.Rn 

Energy-level statistics and wavefunction correlations 
are two basic issues in the theory of disordered electron 
systems at the Anderson metal-insulator transition. So 
far, however, progress here has been largely limited to 
scaling arguments in combination with numerical simu- 
lations - for lack of any parameter at the critical point. 
Our interest in this problem owes much to the recent 
works , where a remarkably simple relation between 
two fundamental characteristics of the critical system, 
the spectral compressibility and the critical exponent of 
the eigenfunction correlations, was derived analytically. 

In the present paper we link the basic idea in [QJ|] 
with the problem of long-term kinetic correlations at 
the metal-insulator transition. Specifically, we address 
the question as to the mechanism of the system-size 
(L) dependence of the ensemble-average critical conduc- 
tance g(L). In accordance with the general principle of 
scale invariance at the critical point, the critical con- 
ductance at the Anderson transition || converges at 
large L to a scale-invariant value g(oo) ~ 1 (in units 
of e 2 /h). What we are interested in is thus the correc- 
tion Sg(L) = g(oo) — g(L), which breaks the macroscopic 
scale invariance and behaves as 



6g(L) ~ (L C /LY 



(1) 



with a critical exponent y > 0. The microscopic length 
scale L c is defined by the condition 5g(L c ) ~ 1 and for 
the conventional Anderson transition in three dimensions 
(d) coincides with the mean free path I. Most of our dis- 
cussion will be focused on the 2d case; namely, on the in- 
teger quantum Hall (QH) critical point. The dissipative 
conductance at the QH transition obeys Eq. (1) with the 
length L c given by ln(L c /Z) ~ g 2 (l); so that, in principle, 
the "elementary block" scale L c may be much larger than 
I. For concreteness, however, we assume in what follows 
that g(l) ~ 1 and L c ~ I. This is the case, e.g., at the 
plateau transition in the lowest Landau level. We seek a 
universal relation between the exponent y and other crit- 
ical exponents, independent of the microscopic details. 
Our main result is the following. We argue that 



L 2 /D 



dtP(0, t) 



(2) 



where P(0, t) is the quantum return probability in time t, 
t the scattering time, D oc L 2 ~ d || the long wavelength 
diffusion coefficient (independent of L in the 2d case). 
The above integral is familiar from the weak-localization 
theory |3) : in the simplest case of time-reversal symme- 
try, the finite-size correction to g 3> 1 in 2d is propor- 
tional to it and grows with L as \n(L/l). At the critical 
point the integral diverges more strongly - in a power- 
law manner. What we show is that at the transition the 
scale-dependent correction to the conductance is given by 
the inverse of the integral. It follows that the exponent y 
is not an independent one, but is equal to the critical ex- 
ponent of eigenfunction correlations rj Q , which controls 
behavior of P(0,t) oc t~ x+r i/ d at large t. 

We start by defining two disorder-averaged two- 
particle correlators 
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M„(r) = — (p e (0,0)p £+w (r,r)) 
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(3) 



where the one-particle spectral function at energy e reads 
p e (r, r') = Y,a S ( € - e c)V'a(i")V'a(i"') and a labels eigen- 
states in a given realization of disorder. Throughout the 
paper we neglect the e dependence of disorder-averaged 
quantities, in particular of the averaged density of states 
p = (p e ), where p e = L~ d J d d rp € (r,r), L d is the sys- 
tem volume. The normalization of the first ( "kinetic" ) 
correlator J d d rP u! (r) — 2ir5(ui) reflects particle number 
conservation, while J d d vM UJ (r) = 27rA _1 (l + R^) is ex- 
pressed in terms of the correlator of spectral fluctuations 
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In the above, A = l/pL d — > is the mean level spacing. 
We now write down the identity 
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where A is an arbitrary parameter characterizing the 
Hamitonian of the system if (A) and oj a p = e a — ep. 
We specify the parameter A by considering d 2 R LJ /dA 2 
a second-order response of the system to the small 
(A — > 0) perturbation SH(A) = Ae lkr . Doing so 
we in fact introduce a k-dependent spectral correla- 
tor R UJ (k,A). Next we define its Fourier transform 
7^ w (r, A) = J ( ^J^ d e^ lkr R LJ (k, A) and expand in powers 
of A. To second order in A [Eq. (5)] the expansion of 
TZu){r) readily yields 

l d 2 MUr) 

A dA 2 7T dcu 2 U 

2 d f duj' 

+ ~p 2 Ad^ J ~ ^ £ (°' r )^+ w '( r '°HP^+"' ~ Pe+^j) • 

We decompose the three-particle correlator in the second 
line of this equation according to the following pattern 

p(dp ei (0,r)S P e 2 (r,0)Spe 3 ) (7) 
- (fyei(0,r)<Spe a (r,0)) ({5p ei 5p €3 ) + (8p £2 6pt 3 )) . 

Note the extreme "non-Gaussian" character of the cor- 
relations: the irreducible three-particle average is rep- 
resented as a product of two pairwise correlators. The 
crucial assumption behind this decoupling is that the cor- 
relations can be separated into fast and slow (in energy 
space) pieces: namely, the spectral correlations [Eq. (4)] 
decay rapidly on an energy scale which vanishes in the 
thermodynamic limit, whereas the eigenfunction corre- 
lations [Eq. (3)] decay slowly according to the diffusion 
law. In effect, if one puts r = 0, Eq. (7) constitutes pre- 
cisely what was incorporated in the Langevin description 
(and was called a "i-space decoupling") in When 
substituted into Eq. (6), it gives, for the last term in the 
r.h.s., 

1 f duj' , . d , 
—r / — P»-{r) — (R u+U1 , - R u - U ,) . 8 
7rA J ur ouj 

The decomposition above is a peculiar property of the 
average (7) and we proceed to discuss the conditions un- 
der which it is legitimate (the conditions will equally ap- 
ply to the method of p],|| ) . If R u would be characterized 
by a single scale on which it falls off sharply, the separa- 
tion of the fast and slow variables would simply require 
that Pu(r) be a smooth function of to on this scale. The 
spectral correlations (in closed systems) indeed start to 
decay on the smallest scale of A; specifically, R w oc oj~ 2 
at A < \u\ < E c @, where E c ~ gA is the Thouless 
energy, and continue to decay at larger |cj|. However, 
the argument that the energy-level correlations are short 
ranged at A — » may be illusive, since R u falls off at 
M si E c in a power-law manner, which necessitates a 
power counting. What really matters to the justification 
of the above decoupling is how fast i? w vanishes at large 
\uj\. E.g., one sees that if R w decays slower than uj~ x , 



all frequencies in the range |u/| < \u\ will give roughly 
the same contribution to the integral over ui' in Eq. (8) - 
in contrast to what the decoupling implies. The sought 
conditions can thus be formulated as follows: (i) R^ falls 
off with increasing |w| not slower than (ii) P^O) 

does not decrease with \uj\ faster than Here and be- 
low we mean the frequency dependence in the diffusive 
regime E c < \w\ < r _1 . 

Clearly, the decomposition in terms of the pairwise cor- 
relators fails in the insulating phase, where is sin- 
gular at u> = [see condition (ii)] even in the thermo- 
dynamic limit. So the real question is about accuracy of 
the decoupling in a system with finite g, where the con- 
dition on P w (r) is not severe but the behavior of R u [see 
condition (i)] becomes crucial. It is worthwhile to em- 
phasize at this point that the validity of the decoupling 
is not related directly to the expansion in powers of 1/g 
and the above conditions can be satisfied (at A/|w| <C 1) 
even if g ~ 1. On the other hand, it is important that 
the diagrammatic analysis 0] up to three-loop order does 
conform to these conditions. Consider particular exam- 
ples. The frequency dependence of R^ becomes steeper 
as the localization effects get stronger and so does it with 
decreasing dimensionality. In the 3d case, at g(l) 3> 1, 
Ru! decays too slowly - as l^l" 1 / 2 |6| - which means that 
the decoupling certainly fails deep in the metallic phase. 
In the critical regime of the metal-insulator transition, 
i? w behaves as \uj\~ x at > A and how x is related 
to other critical exponents is still an open problem [Q. 
However, the remarkable point to notice is that the de- 
coupling which fails in the metallic phase may work at the 
transition, provided x > 1. In 2d, R^ falls off at g 3> 1 
as \lo\~ 1 P|, which is a marginal case in respect of the va- 
lidity of the decoupling (meaning logarithmic accuracy). 
The w _1 behavior of R u at large g implies that x > 1 
at the QH transition, which is one of the above condi- 
tions. Now, we must also be concerned about how (r) 
behaves as u> — > 0. At the transition P u {r) diverges at 
small lj, down to \u\ ~ A, as |cj|~ 1, / rf To decouple the 
spectral and eigenfunction correlations in 2c?, one needs 
rj < 2 (77 = 2 corresponds to an infinitely sparse fractal). 
At the QH critical point ?)~0.4 |). We conclude that 
the decoupling is valid at the QH transition || . 

From now on we specialize to the QH transition. The 
net result of the decoupling (7) is that the expansion of 
7tui{r) i n A for any given u> and r yields just M^(r) and 
Pu>{r). It is convenient to transform to the time repre- 
sentation - introducing M(r,t) and P(r,t) - where the 
relation between the two functions assumes a particularly 
compact form: 

tM(r,t)-2K j\t'P{r,t') = ^p^ . (9) 

Here we have also introduced the dimensionless Fourier 
transform K(t) — A -1 J due~ luJt R U j and a similar trans- 
form IC(r,t) for 1Zu(r). It is worth emphasizing that the 
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validity of the decoupling (the condition (i) above) re- 
quires K(t) to be constant in the thermodynamic limit 
JTo| ]. For our purposes, we thus may write R u = 
K6(u/A) + (A/E)F(u/E), where £ = r" 1 , F(x) is a 
dimensionless function of order unity. It is the last term 
that is crucial to the analysis of the finite-size scaling. 
Clearly, it violates the condition (i); nevertheless, Eq. 
(9) is useful if we are interested only in the critical ex- 
ponents - not the exact (nonuniversal) shape of F(x). 
Indeed, at r = we recover the relation 101 



K= lim -tP(0,t)( [ dt'P(0,t')) 

t/r^oo 2 \Jq J 



(10) 



which, taken at finite t/r, allows to estimate the small 
t-dependent term in K(t). In this way we find that 
K(oo) - K(t) ~ (r/t)''/ 2 (i.e., F(x) diverges as x~ 1+r i^ 
at x — > 0). Note that in going to the last step we 
have neglected the parametric correlations in the r.h.s. of 
Eq. (7) (while in jjj they were omitted by construction), 
which needs explanation. We first rewrite d 2 IC(r,t)/dA 2 
at r = as the response d 2 K(t)/du 2 to the application 
of a white- noise random potential W(r) with the corre- 
lator (H^W^) = u 2 5(r). This additional random potential 
merely renormalizes the scattering rate. It follows that 



d 2 K,(0,t) 
dA 2 



4-7T/9 



dK{t) 
<9£ 



(11) 



(recall that, for simplicity, we keep p independent of S, 
even at g ~ 1). Hence the starting idea in about 
the statistical equivalence of all members of the ensemble 
of Hamiltonians, amounts to ignoring the dependence of 
K(t) on S. One can easily estimate the omitted con- 
tribution to Kit): it is ~ (r/t) 1+n , which is indeed 
smaller than the leading correction obtained above. Now 
we can argue in the same way to make sure that if t 
is large enough, the parametric correlations can be ne- 
glected at finite r as well. The large-i limit means t » 
max{r, r 2 /D}, so that as r increases it needs more time 
for the parametric correlations to be forgotten. It is also 
instructive to see how in 2d the first terms in the expan- 
sion of M, P, and K in powers of g~ Y cancel in Eq. (9). 
Note that M(r, t) and P(r, t) can be split into short and 
long ranged (in real space) parts. Specifically, M(r,t) ~ 
2^/3(0, t) + f 2 (r) P (Q, t) + (27TP)- 1 /„' dt'p(r, f)p(r, t-t'), 
P{r,t) ~ 2npP(r,t) +p(r,t), and K ~ (iirg)- 1 . Here 
(3(r,t) describes ballistic motion at t < t and vanishes 
at larger t, J dt/3(r,t) = / 2 (r), f(r) = e~ r / 21 J (k F r) is 
the Friedel function [kp is the Fermi momentum), and 
p(r, t) is the retarded Green function of the classical dif- 
fusion equation. For concreteness, we mean a system 
with completely broken time-reversal symmetry. Since 
K ~ g~ l , we need next order terms in M(r,t) as com- 
pared to P(r, t) - hence the two-diffuson term above. One 
can see that in the limit t 3> max{r, r 2 /D} all the terms 



cancel out in the l.h.s. of Eq. (9). Despite the simplicity, 
it is a useful check. 

Now we write the correlators in the scaling form which 
they assume at the QH transition: M(r,t) — tjjM(^), 
P{r, t) = 75jP(-j)j), where both M(x) and P(x) behave as 
x^^l 2 at x <C 1 with r\ ~ 0.4 Q. The hydrodynamic form 
of the scaling functions implies that r > ) and ( > t, 
where A is the magnetic length and D ~ p^ 1 ~ A 2 t _1 
(for definitcness we mean the case of short-range disor- 
der in the lowest Landau level) . Clearly, within the scal- 
ing approach the average conductance is scale indepen- 
dent. Indeed, whatever the function P(x) is, the velocity- 
velocity correlator G v (t) = j d 2 rr 2 9 is propor- 
tional, within the scaling description, to 5(t). Hence, 
it is corrections to the scaling that determine the asymp- 
totic behavior of G v (t) at t ^> r. We expect that the 
quantum kinetic correlations decay at the transition in a 
power-law manner: 



(12) 



We now appeal to Eq. (9) in order to find y. The es- 
sential idea is that the tail in G v (t) originates from a 
power-law (in r/t <C 1) correction to the limiting value 
of K(t) at t — > oo. Indeed, at the critical point K(t) 
at A — > is parameterized as a function of the single 
variable t/r. It follows that the r-dependent correlators 
written in the scaling form as functions of r 2 /Dt satisfy 
Eq. (9) only if K (t) is strictly constant. Hence, taking 
into account t-dependent terms in K(t) will generate cor- 
rections to the scaling. As noted above, these yield the 
sought tail in G v (t). We thus expand K(t) in r/t and 
substitute the leading term in the expansion into Eq. (9). 
Multiplying the equation by r 2 and doing the integrals 
over r, we obtain the contribution to G v (t) which is due 
to the dispersion of K(t) - it is given by Eq. (12) with 



y = v 



(13) 



The exponent y defined by Eq. (12) governs behav- 
ior of the average conductance g(L) in a finite sam- 
ple of size L x L. Specifically, the conductance at 
L — > oo is given by the Kubo formula for the conduc- 
tivity 3(00) = 5 p J dtG v (t), and cutting the integral 
over t at t ~ L 2 /D yields g(L). Since we associate the 
long-term tail in G v (t) with the decay of the difference 
K(oo) — K(t), which in turn is given by the behavior of 
P(0,t) [Eq. (10)], we finally arrive at Eq. (2). 

A cautionary remark is in order at this point. We as- 
sumed that the corrections to P(r,t) and M(r,t) do not 
cancel each other in Eq. (9) exactly unless other irrele- 
vant scaling fields are taken into account - since, apart 
from the common exponent y, the corrections depend on 
microscopic details (say the correlation properties of dis- 
order). Eq. (9) thus tells us that there is a contribution 
to G v (i) which scales with y = rj. However, because of 
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the parametric correlations in the r.h.s., we cannot rule 
out solely on the basis of this equation the possibility 
that there exists another contribution with smaller y (a 
"less irrelevant" scaling field), not at all related to the 
behavior of K[t). Therefore, to be precise, what we have 
been able to prove is the inequality y < rj. At this point 
we might seek support from numerical data. These sug- 
gest strongly that the sign < above is in fact "equals": 
indeed, the values of y ~ 0.4 - 0.5 and V - °- 4 

P,|j~H obtained in the numerical simulations are fairly 
close to each other. 

So far we have dealt with the scaling of the average 
conductance g(L). Of what significance could it be in 
real experiments in view of the fact that at the critical 
point the mesoscopic fluctuations are very strong? To 
extract the scale dependent 5g at the QH transition, the 
dephasing length ~ (D/T) 1 / 2 (or, interchangeably, 
Lu, ~ (D/oj) 1 / 2 in ac measurements) should be much 
smaller than the sample size - the measured conduc- 
tance then will be a self-averaging well-defined quantity 
g {Lff,) (independent of L in 2d) with a small tempera- 
ture dependent correction Sg^^L,/,) ~ (X/L^y. In this 
respect, the experiment with the QH metal looks much 
like the measurement of the localization corrections in 
a good metal. Note two points relevant to the experi- 
ment. First, the average ("coherent") conductance g(oo) 
(at L/X — » oo ) and the measured at finite T ("effec- 
tive") conductance g eS (oo) (at both L^/X — > oo and 
L/L^ —> oo) are by no means the same. One way to 
look at the problem is to divide the system at the criti- 
cal point into blocks of size and couple them to each 
other "incoherently" , thus representing the system as a 
classical resistance network which obeys the Ohm's law 
locally. Since the quantum interference is preserved in- 
side the blocks, the elementary conductances gi of this 
network will fluctuate strongly around mean g(L<p). One 
may now consider expanding g cS — g in series in the cu- 
mulants of gi fl3|| - the first term in the expansion is 
given by — ((gj) — (gi) )/2 (gi). Since there is no small 
parameter in this expansion, it follows immediately that 
\g eS (oo) — 3(oo)| ~ 1. What is important to us, how- 
ever, is that the scaling of the corrections to either of 
the quantities is governed by the same exponent y. The 
difference between the measured and averaged conduc- 
tances is nonetheless crucial when referring to the con- 
cept of the universal g at the quantum phase transition 
(in this respect the above arguments complement those 
in [0). Second, the applicability of the non- interacting 
model considered in the paper to the experiment at the 
integer QH critical point requires comment. There exists 
a long-range contribution to the velocity- velocity correla- 
tor G v (t) associated with electron-electron interactions: 
as is well known, at large g in 2d it decays as t^ 1 || - 
i.e., it scales similarly to the weak- localization contribu- 
tion. We have shown that the t^ 1 tail which is due to 



the localization effects is transformed at the integer QH 
transition into i~ 1 ~' ? / 2 . The question as to what extent 
the Coulomb interaction affects this behavior of G v (i) 
requires further study. 

In conclusion, we have studied how the critical conduc- 
tance g(L) at the metal-insulator transition behaves as a 
function of the system size L. We relate the size depen- 
dent correction Sg(L) ~ (l/L) y to the quantum return 
probability and argue that the exponent y is the same as 
•q - the critical exponent of eigenfunction correlations. 
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